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ABSTRACT 

We present Sunyaev-Zel’dovich (SZ) effect measurements from wide-field images towards the galaxy 
cluster RX J1347.5—1145 obtained from the Caltech Submillimeter Observatory with the Multiwave¬ 
length Submillimeter Inductance Camera (MUSIC) at 147, 213, 281, and 337 GHz and with Bolocam 
at 140 GHz. As part of our analysis, we have used higher frequency data from Hersc/iel-SPIRE and 
previously published lower frequency radio data to subtract the signal from the brightest dusty star¬ 
forming galaxies behind RX J1347.5—1145 and from the AGN in RX J1347.5—1145’s BGG. Using 
these five-band SZ effect images, combined with X-ray spectroscopic measurements of the temper¬ 
ature of the intra-cluster medium (ICM) from Chandra, we constrain the IGM optical depth to be 
Te = 7.33to'97 X 10”^ and the IGM line of sight peculiar velocity to be Upec = —lOdOtglo km s“®. 

The errors for both quantities are limited by measurement noise rather than calibration uncertainties 
or astrophysical contamination, and significant improvements are possible with deeper observations. 

Our best-fit velocity is in good agreement with one previously published SZ effect analysis and in 
mild tension with the other, although some or all of that tension may be because that measurement 
samples a much smaller cluster volume. Furthermore, our best-fit optical depth implies a gas mass 
slightly larger than the Chandra-derived value, implying the cluster is elongated along the line of sight. 

Subject headings: galaxies: clusters: intracluster medium — galaxies: clusters: individual: 

(RX J1347.5-1145) 


1. INTRODUCTION 


The peculiar velocities (upec) of galaxy clusters and 
other large-scale structures are related to the underly¬ 
ing cosmology, particularly the properties of dark energy 
and the long length-scale behavior of gravity. As such, 
measurements of Upec can be used to place constraints 


on a variety of cosmological parameters (e.g., Kglsgr 


1987||Scoccimarro|20041|Bhattacharya fc Kosowsky|2Ul)8| 


Kosowsky & Bhattacharya| 2009p. In addition, the 

assumption of large-scale homogeneity in the universe 
based on the Copernican principle can also be tested 
using Upec measurements (Clarkson et al. 2008 Clark¬ 
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son 2012 fPlanck Collaboration et al. 2014). The ki- 
netic Sunyaev-Zel’dovich (SZ) effect, which describes the 
Doppler shift of Cosmic Microwave Background (CMB) 
photons that scatter off of hot electrons in the intra¬ 
cluster medium (ICM), provides an excellent tool for such 


measurements dSunyaev &: Z el’dovich 1972 Birkinshaw 
1999 Carlstrom et al. 2002). In particular, the surface 


brightness of the kinetic SZ effect signal is independent 
of redshift and is directly proportional to the line of sight 
peculiar velocity of the ICM relative to the fixed refer¬ 
ence frame of the CMB. 

However, the relative dimness of the kinetic SZ effect 
signal, along with contamination from several other as¬ 
trophysical sources such as primary CMB fluctuations 
and the thermal SZ effect, present a significant observa¬ 
tional challenge. Recently, progress has been made by 
averaging the kinetic SZ effect signal over large samples 
of objects in wide field CMB surveys. For example, ki¬ 
netic SZ effect measurements based on multi-band data 
from the WMAP and Planck satellites have been used 
to place upper limits on the average and rms peculiar 
velocities of large (:::^ 1000) X-ray selected cluster sam- 


pies (IKashlinskv et al. 12008 

Osborne et al.|2011| |Planck| 

Collaboration et al. 2UI4p. 

In addition, Planck and 


have made separate modest-significance detections of the 
kinetic SZ effect signal using pairwise stack s on galaxies 


detected in the Sloan Digital Sky Surve y (Planck Col- 
laboration et al.||2015a Hand et al.||2012). Furthermore 
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data from Planck, ACT, and the ground-based South 
Pole Telescope (SPT) imply a non-zero kinetic SZ effect 
power spectrum on sr nall angular scales dPlanck Collab- 


instrument description in Sayers et al. 


oration et al.|[2015b[ [George et al.|[^15[ [Sievers et al. 
2013p. Complementary to these large statistical anal- 
yses, several efforts have also been made to detect the 
kinetic SZ effect signal in individual clusters using multi- 


Kitayama et al. [20041 Mauskopf et al.|2012 

Zemcov et al. 

20121 Mroczkowski et al.|2012( [Lindner et a 

.|2015[|Adam 


mg component in the cluster MACS J0717.5-1-37 45 pro- 


et al. 


2013cb 


ducing the h ighest significance detection to date ( [Sayers 
ivsi 


'i'hls analysis focuses on SZ effect measurements to¬ 
wards the cluster RX J1347.5—1145, which is a well- 
studied system that was the target of some of the earliest 


published thermal SZ effect detections (Komatsu et al. 


19991 IPointecouteau et ai.| 

19991 IKomatsu et al.||2UUlt 

Pointecouteau et al.[ 2001 

Reese et al. 2002|). Due to 

its brightness, KX J 1347.5—1145 has also been used tor 
the initial studies published from several recent SZ ef- 

feet instruments (^ason et al. 

2010 

Zemcov et al.[[2012 


SZ effect measurements towards RX J1347.5—1145 have 
been published. 


-1420^1270 km s 


Kitayama et al. (2004) found Cpec = 
based on 21, 150, and 350 GHz 


measurements collected from the Nobeyama 45 m Ra¬ 
dio Telescope and t he Ja mes Glerk Maxwell Telescope, 
and Zemcov et al. (2012) used high spectral resolution 
200-300 GHz data from Z-Spec combined with 140 GHz 
Bolocam data to constrain Cpec = -1-450 ± 810 km s~^p^ 

In this manuscript we present SZ effect measurements 
towards RX J1347.5—1145 collected at 147, 213, 281, 
and 337 GHz using the Multiwavelength Submillimeter 
Inductance Camera (MUSIC) and at 140 GHz with Bolo¬ 
cam. Summaries of the instrument characteristics are 
given in Section]^ and the data reduction procedures are 
described in Section Section details the analysis of 
Chandra X-ray exposures to obtain a spectroscopic con¬ 
straint on the intra-cluster medium (IGM) temperature 
Te. Section discusses the model used to describe the 
shape of the SZ effect signal towards RX J1347.5—1145, 
and Section details how the SZ brightness measure¬ 
ments derived from the model fits, in combination with 
the X-ray spectroscopic measurements, are used to con¬ 
strain the ICM optical depth Te and peculiar velocity 
Cpec. A comparison to previously published results is 
given in Sectionj^ Throughout this manuscript, all phys¬ 
ical quantities have been derived using a fiat cosmology 
with Hm = 0.3, Qx = 0.7, and Hq = 70 km s“^ Mpc“^. 

2. SZ EFFECT INSTRUMENTATION 

2.1. MUSIC 

MUSIG was built to serve as the long-wavelength imag¬ 
ing camera for the Caltech Submillimeter Observatory 
(C SO), and its overall d esign and constr u ction is detailed 
in [Glenn et al. (2008 MaloneyetaL (2010), Golwala 


et al. 


and Sayers et al. (2U14p. Relative to the 


[Kitayama et al.[ ||2004| report a value of -1-1420 km s“^, 
but that value is based on an opposite convention for the velocity 
direct i on co mpared to what is used in this work and in |Zemcov| 
[et al.[|2M^. 


(2014), the cryo¬ 


genic lens has been redesigned in order to change the 
plate scale of the focal plane from 7"/mm to 11.5"/mm. 
As originally built, eight focal plane tiles would fill MU¬ 
SIC’S 14' circular field of view (FOV). However, the ini¬ 
tial fabrication yielded only two working tiles, and this 
plate scale change allows these two tiles to fill a square 
11.5' FOV (14' diagonal). Subsequent to this redesign, 
one of the working tiles was damaged, and all of the data 
described in this manuscript were collected using a single 
tile with an approximate FOV of 5.8' x 11.5' (AZ x EL). 

The MUSIC focal plane uses phased arrays of slot 
antennas to couple to incoming radiation over a broad 


band width spanning from ~ 100-400 GHz (Goldin et al. 
2003). The output from each antenna array is then 
coupled to four separate microwave kinetic inductance 
detectors (MKIDs), with three-pole lumped element fil- 
ters defining a different observin g band for each MKID 
(Kumar et al. 2009[ Duan 2015). These four observing 
bands are centered at 147, 213, 281, and 337 GHz 1^ and 
they have corresponding point spread functions (FSFs) 
that are approximately Gaussian and have full-width half 
maxima (FWHM) of 48", 36", 32", and 29" Each MU¬ 
SIC focal plane tile consists of 72 antenna arrays and 288 
MKIDs, and the tile used to collect the data described 
in this manuscript has 200 MKIDs that pass all of the 
initial quality checks. An additional cut is made based 
on the on-sky sensitivity of each detector, removing all 
MKIDs with a sensitivity more than two standard de¬ 
viations poorer than the median sensitivity within each 
observing band. After this cut 125 MKIDs remain, with 
30, 35, 32, and 28 MKIDs at 147, 213, 281, and 337 GHz. 

2.2. Bolocam 

Bolocam served as the long-wavelength imaging cam¬ 
era at the CSO prior to the installation of MUSIC, and 
the in stru ment is described in detail in IGlenn et al.l 
(1998) and [Haig et al. (2004). Bolocam operated with 
an 8' circular h'UV filled with 144 spiderweb bolome¬ 
ters using neutron-transmutation-doped Ge thermistors, 
of which approximately 115 were optically coupled. The 
data described in this manuscript were collected with 
Bolocam operating at 140 GHz, where the PSFs had a 
FWHM of 58". 

As detailed in Section]^ the Bolocam data provide a 
modest improvement in constraining the SZ effect signal 
compared to using the MUSIC data alone. In addition, 
this is the first SZ analysis using MUSIC data, while the 
Bolocam data for RX J1347.5—11 45 have already been 
published in multip le analyses (e.g., Sayers et al. 2013b|a 


Czakon et al.|2015 ). As a result, the well-vetted Bolocani 


data provide a useful consistency check for the MUSIC 
measurements. 

3. SZ DATA REDUCTION 
3.1. General Reduction and Calibration 

For both MUSIC and Bolocam, the data were collected 
by scanning the CSO in a Lissajous pattern with an rms 

These values represent the effective band centers for a rela- 
tivistically corrected thermal SZ spectrum assuming an ICM tem¬ 
perature of 10 keV. 

These FWHM are slightly larger than the values given in 
[Sayers et al.[ |[2014| due to the optical reconfiguration. 
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TABLE 1 

Bolocam+MUSIC Data 


band (GHz) 

PSF 

HPF (half-signal) 

HPF (2xR250o) 

rms/beam 

image size 

140 GHz 

58" 

(9.5')"^ 

0.89 

1.1 mJy 

14' X 14' 

147 GHz 

48" 

(2.5')-i 

0.30 

1.1 mJy 

12' X 12' 

213 GHz 

36" 

(2.4')“^ 

0.26 

0.8 mJy 

12' X 12' 

281 GHz 

32" 

(2.9+1 

0.33 

2.0 mJy 

12' X 12' 

337 GHz 

29" 

(2.7+1 

0.29 

7.2 mJy 

12' X 12' 


Note. — The Bolocam+MUSIC data used in this analysis. The columns give the 
thermal SZ effect-weighted band center, the PSF FWHM, the angular scale at which 
half of the signal from a cluster-sized object is attenuated due to the effective high- 
pass filtering of the data processing, the fraction of signal remaining for angular scales 
corresponding to 2x_R2500 = 4.1^ (see Section]^, the rms/beam in the image, and the 
total size of the image. 


speed of 240" and an amplitude of 3' and 4' respec¬ 
tively. Nightly observations of Uranus or Neptune were 
made to determine the absolute flux calibration of the 
data, which is estimated to have an rms uncertainty of 
5.0% for Bolocam and 6.0%, 4.6%, 5.5%, and 9.4% in 
the 147, 213, 281, and 337 GHz MUSIC observing bands 
(Sayers et al. 2013a 2014). Included in all of the above 
uncertainties is an overall 3.3% rms uncertainty in the 


normalizat ion of the planetary brightness model (Sayers 
et al.|20l2 ). In the case of Bolocam, hourly observations 


ot bright compact objects were made in order to o btain 
a pointing accuracy of ~ 5" rms (Sayers et al.|2011 ). For 
MUSIC, a single pointing offset derived from the flux 
calibration observations is sufficient to obtain a pointing 
accuracy of ~ 3" rms, and therefore no dedicated obser- 
vations were made for the purpose of pointing calibration 
(jSayers et al.||2014|). 

i'he Bolocam data have been used in several SZ ef¬ 
fect analy ses, and the details of their reduction ar e de- 
scribed in Sayers et al. (2011), Sayers et al.|( 2013a), and 


Czakon et al. (|2U15 ). Briefly, a template for the atmo- 
spheric brightness fluctuations is formed from the aver¬ 
age bolometer signal at each time sample. This template 
is then correlated and regressed from each bolometer’s 
time-ordered data (TOD). Next, a high-pass filter with 
a characteristic frequency of 250 mHz is applied to the 
TOD. Both of these processes attenuate SZ effect sig¬ 
nal, and an image-space transfer function is computed by 
reverse-mapping and processing a cluster model through 
the same reduction pipeline. See Table 
In order to reduce the MUSIC data, the first step 
is to locate discontinuous jumps in the TOD that ap¬ 
pear to be caused by the readout electronics. When 
a discontinuity is found, the TOD for that MKID is 
flagged and removed from the entire 10-minute long ob¬ 
servation. In total, approximately 10% of the data are 
flagged. The MUSIC instrument moves along with the 
CSO telescope during the Lissajous scans, and this mo¬ 
tion relative to the Earth’s magnetic field produces a 
response in the MKIDs at the level of ~ 10 mJy peak- 
to-peak. To remove this response, templates of the AZ 
and EL motion of the telescope are correlated and re¬ 
gressed from the resonant-frequency-response TOD of 
each MKID (changing magnetic fields do not produce 
a signal in the dissipation-response TOD)j^ Additional 


The MUSIC readout works by using a single prob e tone 
centered on the nominal resonant frequency of each MKID ( |Duan| 


noise fluctuations, with 1// and 1//^ power spectra, are 

sourced by gain and phase fluctuations in tne iVIUSIC 

readout electronics. Templates for these fluctuations are 

constructed from electronics fluctuation monitors and 

then subtracted from the MKID TODs. Analogous to 

the Bolocam analysis procedure, a template of the atmo- 

spheric brightness fluctuations is then subtracted, with 

a separate template constructed for each of the four ob- 

serving bands. A high-pass filter with a characteristic 

frequency of 500 mHz is then applied to the I'OD. Com- 

pared to Bolocam, the MUSIC data contain more low- 

frequency noise due to both residual instrumental and 

atmospheric noise, and therefore a more aggressive high- 

pass filter provides an improved signal-to-noise ratio in 

the MUSIC data. Finally, the overall image-space trans- 

fer function of the MUSIC processing is determined sepa- 

rately for each observing band using the same procedure 

as for Bolocam. Due to MUSIC’S slightly smaller F'UV, 

along with the higher frequency high-pass hlter applied 

to the MUSIC TOD, the Bolocam images retain hdelity 

to cluster emission on larger angular scales than the MU- 

SIC images. See Table 1| 






3.2. Point Source Subtraction 


and its emission is non-negligible compared to the SZ 

effect signa l in the Bolocam-tMUSIC observing bands. 

Due to BoTocam-l-MUSIC’s relatively coarse angular res- 


The BCG of RX J1347.5-1145 hosts a strong AGN 


olu tion, and the fact th~ the AON is spatially coi n¬ 

cident with the peak i n the SZ effect signal, it is not 
ossible to model and subtract the AGN emission using 
— ^ ' . As a result, a wide ran^ 


of availaBTe measurements ot the~ACN at frequencies 
close to the Bolocam-|-MUSIC observing bands were used 


to obtain a power-law ht to the hux density of the AGN. 

Specifically, these data incl ude a 1.4 CH z measurement 

from NVSS (ICondon et alTTQQS ), 8.5 an(r22'.5 CHz mea¬ 


sure ments from VLA (Komatsu et al.|2UUl|), a 28.5 CHz 

measurement from (J V K’n7^^^^°lUoE!eeF^m2UU7p 


|et al.|[20Tol l. Changes in the resonance depth and/or frequency 
thus produce a change in the complex-valued transmission of this 
probe tone. A full calibration sweep of the resonance is performed 
once every 30 minutes, and this sweep provides a mapping from 
the complex-valued TOD to separate real-valued TODs that rep¬ 
resent the change in the resonance frequency and resonance depth 
(or dissipation). For MUSIC, the signal-to-noise for astronomical 
sources is many times larger in the frequency response than in the 
dissipation response, and so the dissipation-response TODs are not 
used in this analysis. 
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Fig. 1.— Measurements of the flux density of the AGN 
in the BCG of RX J1347.5—1145 at frequencies close to the 
BolocamH-MUSIC observing bands (see text for more details). 
Overlaid in grey is the 68% confidence region for a power-law fit 
obtained using linmix_err. This fit has a fractional uncertainty 
of 10% at the frequencies corresponding to the Bolocam-hMUSIC 
bands at 140-337 GHz. 


a 30.9 GHz measurement from SZA (Bonamente et al. 


20121), a 90 

GHz measurement trom JVlUS'i'AiNG (jlvo- 

rngut et al. 

20111), 86 and 90 GHz measurements trom 

CARMA fP 

agge et al. 

2013), a 100 GHz measurement 

from NmU 

Komatsu et al. 111999), a 140 GHz measure- 

ment from 

Diablo ( 

Fointecouteau et al. 20011), and a 

250 GHz measurement 

trom yCUHA (Komatsu et al. 


19991. All of these data were collected using mterter- 


ometers except for MUSTANG (which is a high angular- 
resolution single-dish photometer), Diablo (where the SZ 
effect and AGN signals were simultaneously modeled), 
and SGUBA (which was relatively close to the null in 
the thermal SZ effect signal). See Figure H 
fit obtained using linmix.err (Kelly 2007f) 


The best 
which ac¬ 


counts for intrinsic scatter between measurements due 
to potential, unaccounted-for instrumental systematics 
and/or source variability, is 


F{iy) = (59.1 mJy) x 


Extrapolating this fit to the Bolocam-|-MUSlC observing 
bands yields flux densities of 4.37, 4.26, 3.50, 3.03, and 
2.75 mJy at 140, 147, 213, 281, and 337 GHz. Rounded 
to the nearest percent, the rms uncertainty on the extrap¬ 
olation is 10% in all of the Bolocam-|-MUSlG observing 
bands, and this uncertainty is inc luded in the noise esti¬ 
mates as described in Section lOl below. 

In addition to the BCG, the Bolocam-|-MUSIC images 
also contain a non-negligible signal from the population 
of dusty star-forming galaxies behind RX J1347.5—1145, 
although none of these galaxies are individually detected 
at high significance in the Bolocam-|-MUSlC data. The 
emission from these sources has been subtracted accord¬ 
ing to the general procedure de scribed in detail in the 
Appendix of Sayers et al. (2013c I with some slight modifi¬ 
cations described befow. Briefly, a list of candidate galax¬ 
ies is compiled from the three-band Herschel-SPIRE ob¬ 
servations at 600, 850, and 1200 GHz. In total, 192 


TABLE 2 

ffersc/ieLSPIRE Dusty Star-Forming Galaxy Model 


band (GHz) 

140 

147 

213 

281 

337 

total 

normalization 

uncertainty 

-0.73 

1.53 

2.27 

1.58 

1.47 

0.34 

1.12 

0.36 

2.24 

0.80 

1.36 

0.23 


Note. — Bolocam-l-MUSIC-derived fits to the normal¬ 
ization of the dusty star-forming galaxy model obtained 
from Herschel-SPIRE. Averaging all four bands together, 
a best-fit normalization of 1.3(3 it 0.23 is obtained, indicat¬ 
ing that the emission in the Bolocam-|-MUSIC observing 
bands is 36% brighter than the Ifersc/ieZ—SPIRE greybody 
extrapolation. 


galaxy candidates are detected by Herschel-SPIRE with 
a signal-to-noise above 3 in at least one of the Herschel- 
SPIRE bands. For each galaxy candidate, a greybody 
spectrum is fit to the Herschel-SPIRE measurements 
and is extrapolated to the Bolocam-l-MUSlC observing 
bands. For each band, the source positions and bright¬ 
ness estimates for all of the galaxies are combined into a 
single model of the background galaxy emission. Next, 
these models are fit to the Bolocam-|-MUSIC data, float¬ 
ing the value of the overall normalization separately for 
each observing band (see Table [^. The average normal¬ 
ization over the five bands is 1.36±0.23, indicating that 
the background emission from the star-forming galaxies 
is slightly brighter than the greybody extrapolation in 
the Bolocam-1-MUSIC bands. Therefore, the Herschel- 
SPIRE model is multiplied by a factor of 1.36 prior to 
subtracting it from the Bolocam-l-MUSIC data, and the 
rms uncerta inty on this scaling factor is taken to be 0.23 
(see Section 3.3). 

Although a single scaling factor for the Herschel- 
SPIRE model is assumed for all five Bolocam-l-MUSIC 
observing bands, it is possible that the slight mis¬ 
match between the greybody extrapolation and the 
Bolocam-l-MUSIC data is due to the assumed grey¬ 
body shape rather than the overall normalization. How¬ 
ever, the Bolocam-l-MUSIC measurements are not pre¬ 
cise enough to distinguish between the two scenarios, and 
so the simpler case of an overall normalization factor is 
used for this analysis. Furthermore, as described in Sec¬ 
tion the extreme case of not subtracting any model 
of the dusty star-forming galaxies results in only mod¬ 
est changes to the derived ICM constraints, and conse¬ 
quently any minor differences between the assumed and 
true greybody models are unlikely to have a significant 
impact on the overall results. 

Images of the Bolocam-l-MUSIC data, along with the 
background dusty star-forming galaxy and AGN models, 
are shown in Figure 1^ The potential impact of these 
point sources on the SZ effect constraints is explored in 
more detail in Section [6l 


3.3. Noise Estimation 

The Bolocam-l-MUSIC data contain noise from several 
different sources, including the instruments themselves, 
the random arrival of photons, the atmosphere, and the 
astronomical sky. In all cases, the noise is characterized 
by a set of 1000 separate five-band image realizations 
which are used to construct a noise covariance matrix in 
a bootstrap fashion (see Section . Some of the noise 
in each of the 1000 realizations is correlated between the 
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Reduced Data AGN Model SPIRE Model Model Subtracted 



13''47"’36' 30* 24* 13''47'"36* 30* 24* 13''47"’36* 30* 24* 13'’47"’36* 30* 24* 


RA (J2000) RA (J2000) RA (J2000) RA (J2000) 

Fig. 2. — From left to right the columns show thumbnails of the Bolocam+MUSIC data, the model of the AGN emission from the cluster 
BCG, the ffersc/iel-SPIRE-derived model for the brightest dusty star-forming galaxies, and the Bolocam-l-MUSIC data after subtracting 
both of these models. For MUSIC, the solid grey lines show positive S/N in steps of 1, and the dashed grey lines show negative S/N in 
steps of 1 (for Bolocam the lines indicate steps of ±3). The solid black line is at 0. In all cases, the color scale has been stretched to cover 
the maximum and minimum image values. All of the images have been filtered, and the more aggressive high-pass filter applied to the 
MUSIC data results in a more compact size for SZ effect signal compared to the Bolocam data. 


separate bands, and therefore a single noise realization 
contains images for all five bands. Each of these noise 
realizations is constructed as described below, and the 
effects of each noise source on constraining the SZ effect 
signal are detailed in Section 

• Instrument, photon, and atmosphere: A jackknife 
image is formed by multiplying the data in a ran¬ 
domly selected subset of half the observations by 
— 1 prior to producing the image. This jackknife 
preserves the noise properties of the instrument, 
random photon, and atmospheric brightness fluc¬ 
tuations, including correlations between observing 
bands, but it removes all astronomical signals. This 
is the dominant noise in the data, particularly in 
the four MUSIC bands. 

• AGN: The template image of the AGN is different 
for each observing band, and consists of that band’s 
PSF, after accounting for the filtering effects of the 
data processing, multiplied by the estimated flux 


density of the AGN in tha t band. Based on the 
ht presented in Section [3^ a 10% rms uncertainty 
was assigned to the estimated flux density of the 
AGN. Furthermore, because any fractional devia¬ 
tion in the source brightness from the model will 
be common among all five Bolocam+MUSIC ob¬ 
serving bands, the noise fluctuations are assumed 
to be completely correlated between the five ob¬ 
serving bands. Specifically, the AGN template in 
the five bands is multiplied by a single, common 
Gaussian random value drawn from a distribution 
with a mean of 1.0 and a standard deviation of 0.1. 

• ffersc/iel-SPIRE-detected galaxies: The normal¬ 
ization of the 7/ersc/ieZ-SPIRE-derived model was 
determined to be 1.36 ± 0.23, and a Gaussian ran¬ 
dom value based on this constraint is used to de¬ 
scribe the uncertainty on this model. For each ob¬ 
serving band, a template is constructed based on 
the HerschelSVIRE model and that band’s PSF 
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and data processing filter. As with the AGN tem¬ 
plate, the noise fluctuations are assumed to be com¬ 
pletely correlated between all five observing bands 
due to the fact that a fractional change in source 
brightness in one band will result in the same frac¬ 
tional change in the other four bands given the 
nominal greybody fit. 


• Undetected dusty star-forming galaxies: Although 
ffersc/ieZ-SPIRE detects 192 galaxy candidates to¬ 
wards RX J1347.5—1145, the bulk of the popu¬ 
lation of dusty star-forming galaxies is below the 
Hersc/ieZ-SPlRE detection limit. In order to ac¬ 
count for brightness fluctuations due to the non- 
uniform spatial distribution of these galaxies, a 
random sky realization is genera ted based on the 
model of Bethermin et al. (20111, after correcting 
the model tor the galaxies that were detected by 
ff erscfeeUSPIRE and su btracted (see the Appendix 
of ISayers et al.| (|2013cD for more details). A single 
sky realization is used for all five observing bands 
in order to preserve the correlations between bands. 
For each band, the sky realization is convolved with 
the PSF and data processing filter. 


• Primary CMB fluctuations: A random sky realiza¬ 
tion is generate d based on the power spectrum mea- 
surements from Planck Collaboration et al. (2015b I 
and is added to the noise realizations tor each of 
the five observing bands after correcting for each 
band’s PSF and data-processing filter. As with 
the sky realizations of the undetected dusty star¬ 
forming galaxies, a single sky realization is used for 
all five bands in order to preserve the correlations 
between bands. 


• Flux Calibration: The flux calibration uncertainty 
includes two separate components. The overall 
planetary brightness model has an uncertainty of 
3.3%, and fluctuations based on this uncertainty 
are assumed to be completely correlated between 
all five obs ervi ng bands. In addition, as described 
in Section |3.1[ there is an uncorrelated measure¬ 
ment uncertainty on the flux calibration for each 
observing band. 

4. X-RAY DATA REDUCTION 

As described below in Section a Chandra X- 
ray constraint on the ICM temperature Tg is included 
in the overall analysis. This measurement of Tg 
was obtained from the longest Chandra observation of 
RX J1347.5—1145, ID 14407, which was reprocessed us¬ 
ing CIAO version 4.7 (December 2014) and CALDB 4.6.9 
(September 2015). The observation is very clean, with 
62.9 ksec from the 63.2 ksec exposure remaining after 
light curve filtering to reject any 200-second bins where 
the count rate varied by more than 3tT from the me¬ 
dian. The standard ACIS background file was scaled to 
match the off-source count rate in the 10-12 keV energy 
range, and subtracted. Spectral fits were performed with 
the 0.7-9.5 keV band using a photon-absorbed APEC 
plasma model with nn = 4.60 x 10^° cm^ (the Leiden- 
Argentine-Bonn value for RX J1347.5—1145’s sky posi¬ 
tion). Excising the extremely bright cool-core at a ra¬ 
dius of 150 kpc (0.43') results in an emission-weighted 


TABLE 3 

gNFW Model Fit Parameters 


instrument 

band 

(GHz) 

PO 

erg cm“^) 

xVdof 

PTE 

Bolocam 

140 

34.3± 1.5 

147/118 

0.04 

MUSIG 

147 

28.6± 4.0 

239/206 

0.06 

MUSIG 

213 

-16.6±26.9 

456/466 

0.62 

MUSIG 

281 

38.6±11.2 

740/730 

0.39 

MUSIG 

337 

24.8±26.4 

919/971 

0.88 


Note. — Summary of the gNFW model fits to MUSIC and 
Bolocam. T he shape of the mode l is based on the fit published 
in Table 2 of|Czakon et al.| (|2015ll , and the only free parameter 
in the fits is the normalization or the model po. In the limit of 
a purely thermal SZ signal, the value of po will be the same for 
all five observing bands. Except for the 213 GHz band, which is 
near the thermal SZ effect null and easily contaminated by the 
kinetic SZ effect signal, there is good agreement between the 
MUSIC and Bolocam measurements. In addition, the quality 
of fit in all five bands is good, indicating that the data are 
consistent with the assumed model. The PTE values are based 
on the theoretical distribut ion, and therefore m ay be slightly 
biased (see Tables 2 and 3 in|Sayers et al.|2011^. 


Tg = 13.4±0.7 keV within i? 2500 - This measurement 
differs somewhat compared t o previous X-ray meas ure- 
ments of Tg. For example, |Mahdavi et al.||2013 find 
Tg = 12.1±0.4 keV based on a combined Chandra-XMM 
analysis. Furthermore, there is substantial evidence that 


the ICM is not isoth_ermal within i ?2500 (e.g-, Johnson 
et al. 20121. Section explores the impact of these dif¬ 


ferent 1 e values on the overall ICM constraints. 


5. SPATIAL MODEL OF THE SZ EFFECT SIGNAL 


in the ICM to the SE of the cluster core I 

Mason et al. 

2010 Johnson et al.||2012 

Kreisch et al.|20 

H|), its X-ray 

morphology indicates that it is among the 

most relaxed 

known clusters (Mantz et al. 

2015 

). As 

a result, the 

best-tit elliptical generalized INavarro, b’renk, and White 


)2015[), based on a fit to the Bolocam data 


kon et a 

was selected to model the spatial shape of the SZ effect 
signal. This spatial template was then convolved with 
the PSF and data-processing filter specific to each of the 
five observing bands. The normalization of the template 
was varied as a free parameter separately for each ob¬ 
serving band using all of the data within the ove rdensity 


et al.|( 

2015^ 

et al. 

(20101) 


shown in the thumbnails in Figure and the radim pro¬ 
files in Figure Based on these model fits, the average 
SZ effect brightness within R 2500 was computed for each 
band, with the results given in Table 
In performing the model fits, the image pixel weight¬ 
ing and overall values were computed based on the 
assumption that the noise is uncorrelated between pix¬ 
els. Specifically, the image-space pixel sensitivity was 
computed for each observing band based on the standard 


T 2500 is the radius enclosing an average density 2500 times 
the critical density of the universe. As shown in Table [T] the 
Bolocam-I-MUSIG data have good sensitivity to angular scales cor¬ 
responding to R 2500 , allowing for accurate constraints within that 
aperture size. 
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Fig. 3. — From left to right the columns show thumbnails of the Bolocam+MUSIC data after subtracting the A GN and dusty star-fo rming 
galaxy model, the best-fit model of the SZ effect signal based on the elliptical gNFW model from Table 2 of|Czakon et al.| l|2015|l, and 
the residuals after subtracting the SZ effect model. In all cases, the SZ effect model provides a good fit to the data, and the residuals are 
consistent with noise. 


deviation of the pixel flux density weighted by the square 
root of integration time in the set of 1 000 noise realiza¬ 
tions described in detail in Section [3.3| This assumption 
is not strictly correct, as many of the noise terms are 
correlated between image pixels and between observing 
bands. Therefore, to fully account for the noise, the spa¬ 
tial template of the SZ effect signal was also fit to each 
of the 1000 noise realizations and these fits were used 
to compute the full 5x5 covariance matrix for the SZ 
effect brightness. The covariance matrix is given in Ta¬ 
ble 0 and the relativ e co ntributions of the noise terms 
described in Section 13.31 to this covariance matrix are 
plotted in Figure Many of the off-diagonal elements 
of the SZ effect brightness covariance matrix are com¬ 
parable in magnitude to the on-diagonal elements, par¬ 
ticularly for adjacent observing bands. In the case of 
MUSIC, which obtains simultaneous multi-band images, 
some of this correlation is due to fluctuations in the at¬ 
mospheric brightness. However, the bulk of the band-to- 
band correlations are due to the unwanted astrophysical 
signals and the common planetary calibration model. In 


addition, primary CMB fluctuations are not a significant 
noise term in any of the four MUSIC observing bands due 
to the high-pass filter applied when processing the data, 
which corresponds to a characteristic angular multipole 
of 4 - 7000-9000 (see Table 0. 

A range of other models from the literature were also fit 
to the data in order to explore how the SZ effect bright¬ 
ness might change as a result of the chosen spatial tem- 


the average profiles from 

Arnaud et al. (I20I0I, 

Planck 

Collaboration et al. 

(201, 

II, and 

Sayers et al. 

(20I3a| 

along with the KX J 1347.5—1145 specihc proh 

le from 


Bonamente et al. (20121. None of these models produced 
an acceptable ht quality to the 140 GHz Bolocam data, 
with probabilities to exceed (PTEs) below 0.001 in all 
casesj^ As a result, none of these models are suitable 


18 -pjjg pQQ]. fit quality of the ensemble-average p rofil es from 
Arnaud et al. (|2010ll,|Planck Collaboration et al.|||2013ll, and|Sayers| 
et al.| l|2Ui3a ' to the Bolocam data is not surprising, given the 
large amount of cluster-to-cluster scatter relative to those average 
profiles. In contrast, and analogous to the nominal [Czakon et al.| 
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Fig. 4.— Radial profiles of the SZ data and the elliptical gNFW model used to describe it. In order to compare the model to the 
data, it must first be convolved with the transfer function of the data processing to determine how it will appear after this filtering. The 
best-fit model to the data is shown in black, which a solid line representing the model prior to filtering and the dot-dashed line representing 
the model after filtering. The filtering attenuates signal on large angular scales, resulting in a peak with smaller absolute brightness. In 
addition, the filtered signal passes through zero at intermediate radii, and then takes on the opposite sign at large radii. The data, with 
error bars, are shown in red and produce a good fit quality to the model in all five bands. The location of R 2500 is shown as a blue dashed 
line. 


TABLE 4 

SZ Effect Brightness 


band (GHz) 

140 

147 

213 

281 

337 

MJy sr“^ X 10“^ 

-17.7 

-14.3 

2.9 

18.2 

18.2 

Covariance [MJy sr ^ 

|2 X 10' 

-4 


140 GHz 

1.3 

0.3 

0.1 

-0.3 

0.0 

147 GHz 

— 

4.8 

0.4 

-0.3 

-1.5 

213 GHz 

— 

— 

21.7 

3.6 

7.3 

281 GHz 

— 

— 

— 

27.9 

7.0 

337 GHz 

— 

— 

— 

— 

374.4 


Note. — The average surface brightness within R 2500 
derived from the SZ effect model fits. The top row gives 
the best fit surface brightness and the bottom set of 5 x 5 
entries gives the covariance matrix derived from fits to the 
1000 noise realizations. 

for exploring the model dependence of the results. 


( |2015[ l profile used in this analysis, the iBonamente et al.l (|2012|| 
pro tile is based on a fi t speci fically to RX J1^47.5—1145. However, 
the IBonamente et al.l' 2012|) profile is more centrally peaked than 
thefuzakon et aiT|(|2Ulo 1 protile, an d this difference is the mos t likely 
cause ot the poor lit quality of the|Bonamente et al.|(|2012|| profile 
to the Bolocam d ata. The reason to r th e ditterence between the 
profiles derived by |Czakon et al.] | |2015[ | an d|Bonamente et al.|(|2012|| 
is unclear, althougJi several possibilities exist, sucJi as ditterences 
in measurement technique, observing frequency, and point-source 
treatment. 


D. SZ, i i T l i S 

The total SZ effect brightness can be described accord¬ 
ing to 

I{li) =/(i/)tg2 + /(,^)kSZ 
I{^)tSZ = f{i',Te)h{iy)y 

I (j^)kSZ = -RkSZ (iz, Te)h(iy)^^Te 

c 

ksO'T 


y = 


rn. 


^Tfdl 


niecz 
Te= (JTnedl 


where v is the observing frequency, ks is Boltzmann’s 
constant, gt is the Thompson cross section, me is the 
electron mass, c is the speed of light, Ue and Te are the 
electron density and temperature, fpec is the ICM line 
of sight peculiar velocity, f{u,Tf,) describes the spectral 
dependence of the thermal SZ effect signal including rela¬ 
tivistic corrections, h(iz) describes the conversion from a 
CMB temperature fluctuation to a brightness fluctuation 
in MJy sr“^, i?ksz(i^, T?) describes the relativistic correc¬ 
tions to the kinetic SZ effect signal, and dl is the differ¬ 
ential element along the line of sight. Based on these 
equations, the total SZ effect signal is dictated by four 
physical properties of the ICM; Tg, y, Te, and Upec- Al¬ 
though the five-band SZ effect brightness measurements 
described in Sectionj^could in principle constrain all four 
of these ICM properties, a prior for the value of Tg is em¬ 
ployed based on the X-ray measurements from Chandra 
described in Section]^ In addition, the ICM is assumed 
to be isothermal, which implies that y oc TgTe. 

A Markov chain Monte Carlo (MCMC) is employed to 
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Fig. 5. — The total covariance matrix is the sum of several noise terms, as described in Section [3.3| To illustrate their relative contribution 
to the total noise, the plots above show the covariance matrix of each noise term normalized by the diagonal elements of the total covariance 
matrix. Specifically, where is the total covariance matrix and is the covariance matrix of noise term x. The plots 

above show Cf 1 with the total noise covariance (top left), the Herschel—SPTRE dusty star-forming galaxy template (top right), 

noise in the jackknife realizations that is not correlated between image pixels (middle left, approximately equal to the raw photon and 
instrumental noise), the additional noise in the jackknife realizations that is correlated between image pixels (middle, approximately equal 
to the residual atmospheric noise), flux calibration uncertainties (middle right), dusty star-forming galaxies below the iTersc/ieZ—SPIRE 
detection limit (bottom left), uncertainties in the model for the emission from the central AGN (bottom middle), and primary CMB 
fluctuations (bottom right). Note the different color scales for the three rows. For MUSIC, the effective high-pass filtering of the data 
processing has a characteristic scale of £o ~ 8000, and this filter eliminates nearly all of the primary CMB fluctuations. 


sample the parameter space in and ttpec allowed by 
the five-band SZ effect brightness measurements given in 
Table 1^ along with the X-ray prior on Tg. For each step 


in the chain, the SZ effect spectrum is com 

puted using 

the SZpack code described in 

Chluba et al. 

(2012 

20131 

with the '3D' setting, 'i’he results ot the IVICIVIC are 


presented in Figure!^ which shows the clear degeneracy 
between the values or Tg and Upec- Separately marginal¬ 
izing over each parameter yields 68% confidence regions 
of Tg = T.SS^q;^® X 10“^ and Upec = — lOdOtg^o km s“^. 
A plot of the 68% confidence regions for the thermal, ki¬ 
netic, and total SZ effect brightness is shown in Figurej^ 
with the five-band brightness measurements overlaid. 

In order to test the consistency of the SZ effect bright¬ 
ness measurements in the five separate observing bands, 
along with providing an understanding of the contribu¬ 
tion of each observing band to the overall constraints on 
Tg and Upec, five additional fits were performed after dis¬ 
carding the data from one of the observing bands. In all 


cases, the best-fit values of Tg and Wpec shift by less than 
the 68% uncertainties, and in most cases the shift is much 
less than the uncertainty value. In addition, the uncer¬ 
tainties on Tg and Upec generally do not increase signifi¬ 
cantly when one band is dropped, although dropping the 
281 GHz data results in uncertainties that are approxi¬ 
mately 50% larger for both Tg and Upec- This degradation 
occurs because the 140-213 GHz measurements do not 
provide a sufficient spectral range to fully disentangle the 
thermal and kinetic SZ effect signals, and the 337 GHz 
data exhibit significant noise and do not provide a strong 
constraint. In addition, removing the 140 GHz Bolocam 
data, which are approximately two times more sensitive 
to the SZ effect than the 147 GHz MUSIG data, only 
results in a ~ 15% degradation to the Upec constraint. 
This is due to the fact that the separation of the thermal 
and kinetic SZ effect signals is limited almost entirely by 
the uncertainty on the 281 GHz measurement. 

To estimate the systematic uncertainty associated with 
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Fig. 6.— Posterior distributions from the MCMC for the optical depth r© and line of sight velocity Upec derived from the five-band SZ 
effect brightness measurements. Light green indicates the 95% confidence region and dark green indicates the 68% confidence region. The 
marginalized 68% confidence regions are r© = 7.33itQ'gy x 10“^ and Upec = —1040i*^g^Q km/s. 



frequency (GHz) 

Fig. 7.— Five-band SZ effect brightness measurements within 
7^2500 towards RX J1347.5—1145. The error bars denote the square 
root of the diagonal elements of the covariance matrix. Overlaid 
in orange, green, and blue are the 68% confidence regions for the 
thermal, kinetic, and total SZ effect brightness. 


TABLE 5 

ICM Optical Depth and Velocity Constraints 


140 

147 

213 

281 

337 

Te 

^pec 

GHz 

GHz 

GHz 

GHz 

GHz 

10-3 

km s~^ 

/ 

/ 

/ 

/ 

/ 

•7 oo"l"0.96 

• "^’^-0.97 

1 

0 

0 

1 + 

00 00 

0 0 


/ 

/ 

/ 

/ 

’^•^^-1.13 

-13001““ 

/ 


/ 

/ 

/ 

7 Q 4 +O .95 
G-34_o.95 

-890l 

/ 

/ 


/ 

/ 


920 

/ 

/ 

/ 


/ 

7.87tlTr 

i/ciU_ii 8 o 

/ 

/ 

/ 

/ 



— 1130^ 940 

ilOU_ g 3 Q 


Note. — Marginalized 68% confidence regions for r© and 
Up©© based on fits to the full data set (top row) and to subsets 
of the data. The check marks indicate which observing bands 
were included in the fit. Removing the data from one observ¬ 
ing band does not cause the best-fit values to move outside of 
the 68% confidence regions, although the uncertainties become 
significantly larger when the 281 GHz band is excluded. Due 
to the relatively high noise in the 337 GHz band, removing the 
281 GHz band significantly reduces the spectral coverage of the 
data and makes it more difficult to separate the thermal and 
kinetic SZ effect signals. 


the X-ray derived value of for RX J1347.5—1145, fits 
were also performed based on different X-ray analyses. 
First, if the core region is not excised from the Chan¬ 
dra data, then the best-fit value of Te is 12.7 ± 0.3 keV 
(compared to 13.4 ± 0.7 keV when the core is excised). 
This value of Te results in Te = 7.52tlo6 ^ 
fpec = —740^970 km s“^, which differ by 0.19 x 10“^ 
and 300 km s“^ from the measurements based on a 
core-excised T^. An additional fit was also p erformed 
using the va lue of Te = 12.1 ± 0.4 keV from 


Mahdavi 


et al. (20131, which is based on a joint Chandra/XMM 
analysis. 'Ihe best fit values of Te and Upec using the 


Mahdavi et al. (2013) measurement are 7.89lj g 2 ^ 10“^ 
anT"^770^|^km s“^, which differ by 0.56 x 10“^ and 
270 km s“^ from the values obtained using the nomi¬ 
nal measurement of Te. Compared to the measurement 
uncertainties, the differences due to the X-ray analysis 
are smaller by a factor of approximately two for Te and 
by a factor of approximately three for Upec, indicating 
that the range of values for Te obtained from different 
analyses does not significantly change the derived ICM 
constraints, particularly for Upec. 

In addition, to determine the effect of assuming an 
isothermal ICM, rather than a multi-temperature ICM, 
a fit was performed base d on the measurements presented 


Johnson et al. (2012). Specifically, they find two hot 


regions within the ICiVf with Tg ~ 25 ± 5 keV. Based on 


the high-res olution SZ effect results presented in Mason 


et al. (2010), these hotter regions are likely to contribute 


10% of the total SZ effect signal. If 90% of the ICM is 
assumed to have Tg = 13.4±0.7 keV, and 10% is assumed 
to have Tg = 25 ±5 keV, then the best fit values of Tg and 
Vpec are 6.55toJ5 x 10“^ and —1140^900° km s“^, which 
differ by 0.78 x 10“^ and 100 km s“^ from the values 
obtained in the nominal isothermal fit. Particularly for 
Vpec, these differences are smaller than the measurement 
uncertainties, and indicate that the isothermal ICM as¬ 
sumption does not have a strong effect on the derived Tg 
and Vpec constraints. 

As described in Section |3.2[ a wide range of lower 
frequency observations were used to model the emis¬ 
sion from the central AGN and to subtract it from the 
Bolocam-|-MUSIC data. To estimate the potential im¬ 
pact of the AGN emission on the final results, a fit 
was performed without the AGN subtracted from the 
Bolocam-l-MUSIC images. The associated best-fit values 
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of Te and Vpec are 7.39^^ 02 ^ —1910^840 km s“^, 

which differ from the nominal best-ht values by 0.06 x 
10“^ and 870 km s“^. In the case of Upec, this differ¬ 
ence is significant, and indicates that the potential bias 
from the AGN emission is comparable to the measure¬ 
ment uncertainties. Analogously, a Ht was performed 
without subtracting the 77ersc/ie/-SPIRE-derived model 
of the dusty star-forming galaxies, and the best-fit val¬ 
ues of Te and Upec from this fit are 6.9lti;o4 x 10“^ and 
—560^ggg km s“^. These differ from the nominal best- 
fit values by 0.42 x 10“^ and 480 km s“^, and indicate 
that the potential Upec bias from the dusty star-forming 
galaxy emission is approximately 2 times lower than the 
potential bias from the AGN emission while the potential 
Te bias is larger. 


7. COMPARISON TO PREVIOUS RESULTS 

Our derived Upec = — lOdO^g™ km s ^ is in good 
agreement with the previously published kin etic SZ ef- 
fect value of Upec = km s“^ from 


Kitayama 


et al. (20041. Furthermore, both of these measurements 


are consistent with the expected distribution of the bulk 
component of Upec for massive clusters, which has an 
rms of ~ 2 00-30 0 km s“^(e.g., Hernandez-Monteagudo 


& Sunyaev 2010). However, both results are noticeably 


offset from the value of Upec = -1- 450 ±810 km s ^ de - 
rived in the kinetic SZ analysis of Zemcov et al. (2012), 


although there are a number of systematic dinerences 
in measur ement technique that could explain this differ¬ 


ence. The Zemcov et al. 


(2012) measurement was made 
using the single-beam Z-Spec instrument and was there¬ 
fore only sensitive to the value of Upec within a ~ 30" 
diameter region near the cluster center (comp ared to the 
more extended bulk cluster measurements in |Kitayaina] 


et al. (2004) and in this work). In relatively small regions 
of the IGM, similar to what was probed in the |Zemcov| 


et al. (2012) analysis, the inter nal IGM motions can be 
significant (cs 500 km s“^, e.g., Nagai et al.||2003 Dolag 


& Sunyaev 2013), and such motions could easily recon- 
cife the results. In particular, merger activity near the 
center of RX J1347. 5 —1145, for which evide nce is seen by 


Mason et al . (2010), Johnson et al. (2012), and Kreisch 


et al. 12014), could potentially source such large internal 
motions m that region. 


Our derived Tg = 7.39lj;o2 x 10“^, which gives 
the mean line-of-sight integrated electron optical depth 
within the plane-of-sky aperture i? 2500 i can be used 
to obtain a constraint on the gas mass (Afgag) of 
RX J1347.5—1145. First, the parametric model given 
in Section]^ is used to convert between Tg and Tg^sphere = 
5.48 ± 0.77 X 10-3, where Te,sphere represents the mean 
optical depth integrated within the line-of-sight defined 
by 7?250 o and the small asymmetry in the uncertainties 
on Tg is neglected. The gas mass can then be determined 
according to 


Afgas — 


^e,sphere \ 


(JT 




(t-^250o) (MeTOp) , 


where /Xg = 1.14 is the mean molecular weight per free 
electron and rUp is the proton mass, yielding Afgas = 
11.9 ± 1.7 X 10^3 compare our SZ-derived Mgas 

to prior X-ray-derived Mgas values in Table including 


the value of Mgas found by Gzakon et al. (2015) as part 


of their X-ray analysis to detme the value ot R 2500 used 
throughout this manuscript. In each case, the SZ estima¬ 
tion of Mgas was made using the R 2500 value obtained by 
the indicated X-ray analysis. The prior on Tg was kept 
at its nominal value, as small changes in aperture have a 
negligible effect on Tg relative to its uncertainty. 

This comparison indicates that, on average, the anal¬ 
yses incorporating SZ data yield systematically higher 
values of Mgas than the X-ray-only analyses. Such a dif¬ 
ference could easily be explained by an elongation of the 
cluster along the line of sight, which would boost the 
SZ signal co mpared to the X-ray signal (e.g., |Morandil 
et al. 2012). This result is i n good agreement with 


Pointecouteau et al. (2001) and Kitayama et al. 
who also hnd a slight preference for line-of-sig. 


(2004) 
it elon¬ 


gation based on an SZ/X - ray comparison. I n con trast, 
Chakrabarty e t al. (2008), Bonamente et al. (2012), and 
Plagge et al. 12013) all find signihcantly lower SZ sig- 
nals for RX J 1347.5—1145 than expected from the X-ray 
data, and those analyses imply a line-of-sight compres¬ 
sion with axial ratios ranging from ~ 3-5. The reason 
for the disagreement between these two sets of results is 
not clear, although it m ay be related to differ e nces i n the 


SZ data. For exam ple, Pointecouteau et al. (2001), Ki¬ 
tayama et al. (2004), and this work use bolometric imag¬ 


ing da t a at mm wavelengths w hile 
(2008), Bonamente et al. (2012), and 


Ghakrabarty et al. 

Plagge et al. 


use interferometric data at cm wavelengths. In addition 
all six SZ datasets have sensitivities to different angular 
scales. Furthermore, the techniques used to ht the SZ 
and X-ray data, and to compare them, are different in 
all six of the analyses. Resolving these discrepancies is 
beyond the scope of this work and will likely require a 
dedicated, uniform re-analysis of all the relevant data. 


8 . SUMMARY 

By combining four-band MUSIC images with a 
140 GHz Bolocam image we have constrained the thermal 
and kinetic SZ effect signals towards RX J1347.5—1145. 
Based on these constraints, and including a prior 
on the IGM temperature derived from X-ray spectro¬ 
scopic measurements from Chandra, we hnd Upec = 
-1040t3™ km s 3. This value is in good agreement 
with t he previously published result from Kitayama et al. 

( 2004|) , but is i n mild tension with the velocity found by 
Zemcov et ai. (2012). However, some or all of that ten- 
sion may be due to internal cluster motions within the 
relativ ely small cluster volume sampled by [Zemcov et al!) 
(2012). We hnd a mean electron optical depth within 
R2500 of Te = 7.39ti;o2 X 10“3, which implies a gas mass 
within i?2500 of Mgas = 11-9 ± 1.7 x 10^3 Mq. This 
mass is slightly higher than previous values of Mgas de¬ 
rived from X-ray analyses, implying a slight elongation 
of RX J1347.5—1145 along the line of sight. 

For this analysis, we have included a rigorous treat¬ 
ment of the photon and instrumental noise, noise due to 
huctuations in the atmospheric brightness, uncertainties 
in the planetary hux calibration model, and noise from 
three unwanted astrophysical sources (primary CMB 
huctuations, dusty star-forming galaxies, and the clus¬ 
ter’s central AGN). In the case of the dusty star-forming 
galaxies, we have used data from Herschel-SPIRE to sub- 
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TABLE 6 

SZ AND X-RAY Derived Mgas 


X-ray analysis 


^2500 (Mpc 

X-ray Mgas (10^® Mq) 

SZ Mgas (1013 Mq) 


Czakon et al.l 

IMs 

■ 

0.71 

9.2 ± 1.0 

11.9± 1.7 

1 

uonanue et ai. 


1 

0.74 

9.5 ±0.5 

11.7± 1.8 

iviantz et al. ( 

2DTir 

0.80 

13.6 ± 1.5 

13.1 ± 1.8 


Note. — Comparison of Mgas values obtained from three prior X-ray-only analyses 
to that obtained from this analysis incorporating SZ data. For each comparison, the 
SZ analysis is redone using the indicated X-ray-derived value of i?2500i though with a 
common Te prior as described in Section!^ As noted in the text, variations in Te due to 
these small changes i n aperture are neeligl ble c ompared to the measur ement uncertainty 
on Tg. For both the|Czakon et al.|||2015|l and|Donahue et al.|||2014|l comparisons, the 
SZ data find a value ot iMgas that is appr oximately '2a l arger than the X-ray value, while 
the agreement is good in the case of the|Mantz et al^||2014|l comparison. 


tract the brightest objects. Most of the noise sources 
are correlated between the five observing bands, and 
all of these covariances have been accounted for in our 
overall results. For the four MUSIC bands, the total 
noise is largely due to the combination of photon, in¬ 
strumental, and atmospheric noise, which implies that 
longer integration times would yield significantly better 
constraints. For example, if we eliminate these noise 
sources from our analysis, then the uncertainties on Te 
and Upec are reduced by approximately a factor of two to 
±0.44 X 10“^ and ±460 km s“^, respectively. In addition, 
more advanced data-processing algorithms to remove the 
atmospheric fluctuations that are correlated between the 
four MUSIC observing bands may produce noticeab le im¬ 
provements (Adam et al.||20I4 Benson et al.||20d3). 
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